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ABSTRACT 

In this paper we consider the problem of using reduced order dynamic com- 
pensators to control a class of nonlinear parabolic distributed parameter sys- 
tems. We concentrate on a system with unbounded input and output operators 
governed by Burgers’ equation. We use a linearized model to compute low- 
order-finite-dimensional control laws by minimizing certain energy function- 
als. We then apply these laws to the nonlinear model. Standard approaches 
to this problem employ model/controller reduction techniques in conjunction 
with LQG theory. The approach used here is based on the finite-dimensional 
Bernstein/Hyland optimal projection theory which yields a fixed-finite-order 
controller. 
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1 Introduction 


In recent years considerable attention has been devoted to the problem 
of using feedback to control fluid dynamic systems. This problem is com- 
plex and particularly difficult when one is faced with phenomena such as 
shocks. Moreover, these systems are governed by nonlinear partial differen- 
tial equations so that the natural state of the system is infinite dimensional. 
If one assumes that “full state feedback” is necessary to design practical 
controllers, then one would conclude that feedback control of fluid dynamic 
system is “not practical” . However, it is well known that even in finite di- 
mensional control systems one rarely has the ability to accurately sense all 
states, so that some form of dynamic compensation must be used. 

This idea clearly extends to infinite dimensional problems and there is 
a growing literature on observers/compensators for distributed parameter 
systems. In this paper we consider a boundary control problem governed by 
Burgers’ equation. We selected this problem because Burgers’ equation is 
an infinite dimensional model that captures some phenomena (e.g., shocks) 
often observed in fluid flows and because it is simple enough to provide 
real insight into the problem. The goal is to show that it is possible to 
use modern control theory to produce practical finite dimensional dynamic 
compensators for boundary control of nonlinear partial differential equations 
of the type that occur naturally in fluid dynamics. 

We shall present a short summary of one approach (the optimal projec- 
tion method due to Bernstein and Hyland) and show how this approach can 
be used in conjunction with standard numerical schemes to produce a real- 
izable low order controller. The optimal projection method is one of many 
approaches to this problem. However, we shall concentrate on this method 
because a very nice theory has already been developed (for bounded input 
and output operators) and we are more interested in illustrating (to non- 
experts) that recent results in distributed parameter control theory can be 
used to design practical feedback laws, than in discussing the “best” ap- 
proach to the problem. It will be clear from our presentation that we are 
writing for those that are not necessarely “control experts” . The extension 
of the general theoretical results to unbounded input and output operators 
will appear in a forthcoming paper. However, for the compensators pre- 
sented here, we do not need the most general theory since we use the finite 
dimensional version of the optimal projection method. 

As noted above it is almost impossible to observe the whole state. Con- 



Optimal Fixed-Finite-Dimensional Compensator 


2 


trols and sensors are limited to a few points or segments of the boundary, so 
it is necessary to construct an appropriate observer (estimator) of the state 
and design a feedback control law (called a compensator) based on the in- 
formation available from the observed (estimated) state variable. Boundary 
control and observation often leads to unbounded input and output opera- 
tors. Stabilization by dynamic feedback or compensation has been consid- 
ered by Curtain [5], Fujii [8], and Nambu [12] for classes of parabolic as well 
as hyperbolic systems, including control and observation at the boundary. 
All of these approaches produce stabilization schemes that either have the 
same finite order as that of a high-order approximate model, or alternatively, 
open-loop model reduction or closed-loop control reduction techniques are 
applied to achieve a lower-order compensator. An advance was made by 
Schumacher [15], when he gave a theory for designing finite-dimensional 
compensators for a large class of systems, including parabolic and delay sys- 
tems. However, in his theory it was assumed that the control and observation 
operators are bounded. Curtain [4] presented an alternative compensator 
design which applied to the same class of systems, except that unbounded 
inputs and outputs were allowed. In [14], Pritchard and Salamon established 
a framework based on semigroup theory for treating the linear quadratic reg- 
ulator problem for infinite-dimensional systems with unbounded input and 
output operators. Their approach is based on a weak formulation of the Ric- 
cati equations which characterize the optimal feedback law in an appropriate 
dual space. 

Here we consider the problem of designing a fixed-finite-dimensional com- 
pensator for a class of distributed system governed by Burgers’ equation, 
where the control find the observation are implemented at the boundary of 
the domain. The possibility of applying this approach to distributed pa- 
rameter systems was first suggested by Johnson in [9] and Pearson [13]. 
The idea of fixing the order of the finite-dimensional compensator, while 
retaining the distributed parameter model was expanded and developed by 
Bernstein and Hyland in [1] and [2]. The method extends the full order 
LQG case to an “optimal fixed-finite-order compensator” characterized by 
four equations; two modified Riccati equations and two modified Lyapunov 
equations, coupled by an oblique projection whose rank is precisely equal to 
the order of the compensator. Bernstein and Hyland assumed that the con- 
trol and observation operators were bounded and hence boundary control 
and observations were not covered by their theory. 

We will present a Bemstein/Hyland type fixed-finite-dimensional com- 
pensator design, which does extend to unbounded input/output problems. 
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In Section 2 we discuss the existence of a finite-dimensional compensator for 
parabolic distributed parameter systems with unbounded control and obser- 
vation. In Section 3 we summarize the infinite-dimensional optimal projec- 
tion theory from [1], and derive the corresponding equations and feedback 
gains which characterize the fixed-finite-order compensator. In Section 4 we 
present an example, construct the approximation schemes, and discuss the 
computational algorithm used for the optimal projection design synthesis. 
Finally, Section 5 contains numerical results and Section 6 is devoted to a 
few closing remarks. 

2 A Theoretical Existence Result 

We consider the following abstract Cauchy problem 

z(t) = Az(t) + Bu(t), z(0) = z Q G H (1) 

y(t) = Cz(t) t> 0 (2) 

where H is a Hilbert space, u(.) G L 2 ( 0 T;5R m ), y(.) G L 2 (0 T;5ft*), and 
A is the infinitesimal generator of analytic semigroup S(t) on H, generally 
unstable, with exponential growth rate 

W ° = log > 0 ( 3 ) 

so that 

||S(*)|| £(// ) < Me (wo+i ^ for all e > 0, t > 0 (4) 

for some constant M = M(wq,c) > 1. Throughout the remainder of this 
paper we let A denote the translation A = —A + wl, where w is fixed 
and w > wo, so that A has well-defined fractional powers (A)** on H and 
— A is the generator of a strongly continuous analytic semigroup S(t) on H 
satisfying 

\\S(t)\\ C (H) < Me'*', t>0. (5) 

In order to allow for unbounded operators B and C, we assume that B 
G V) and C G C(W , 5ft*), where W and V are also Hilbert spaces such 

that 

V(A) CW^H^V (6) 


with continuous dense injections. More precisely, we assume that B * is 
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[A*] 7 -bounded, or equivalently, 

[A]~^B 6 £(» m , H ) for 0 < 7 < 1. (7) 

Similarly, for the operator C we assume that 

C[A}-' 1 € C(H, 3^) for 0 < 7 < 1. (8) 

It is helpful to interpret (1-2) in mild form. In particular, the solution 
z(t) is given by 

z{t) = S(t)z 0 + [ S(t — s)Bu(s)ds 1 0 < t < T (9) 

Jo 

and the output by 

y(t) = CS(t)zQ + C f S(t — s)Bu(s)ds. (10) 

Jo 

We assume that S(t) is also an analytic semigroup on W and that the 
following hypotheses are satisfied: 

(H-l) There exists a constant b(T ) > 0 such that for every T > 0, 
/ 0 T S(T - s)Bu(s)ds 6 W and 

|jT r 5(T- S )B«( S )ds||^ < 6(r)||u(-)ii L2(0 (11) 

for every u(-) € L 2 (0 T; 3f? m ). 

(H-2) There exists a constant c(T) > 0 such that for every T > 0, 

J l|C'5(t)i|| £ i (0 T; #t)dt < c(T , )||a:|| v (12) 

for every x G W. 

We now give sufficient conditions which imply that the system (1-2) can 
be stabilized by a finite-dimensional compensator of the form 

u/(0) = w 0 


w{t) 

u(t) 


A c w(t) - B c y{t) 
C c w(t) 


(13) 

(14) 
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where A c E R Ne x , B c E 3l NcXTn , and C c E are suitably chosen 

matrices. We need the following well-posedness result for the connected 
system (1-2) and (13-14). This result and proof may be found in [4]. 

Proposition 2.1 Let (H-l)-(H-2) be satisfied, then for all z Q E W, w 0 E 
R Nc there exists a unique solution pair z(t) and w(t) of (1-2) and (13-14)- 
This means that z(t) is continuous in H and absolutely continuous in V, 
that (1) is satisfied for almost every t > 0 where u(t) is given by ( 14 ), and 
that w(t) E R Ne is continuously differentiable and satisfies (13) where y(t) 
is given by (2). 

In addition to hypotheses (H-l) and (H-2), we assume: 

(H~3) Stabilizability Condition (S.C.) 

There exists an operator F G C(H, 3? m ) such that Ap = A + BF 
generates an analytic semigroup Sp(t) = and Sp(t) is expo- 

nentially stable on H, i.e., 

\\SF(t)\\ C[H) < Mpe for wp > 0. (15) 

(H-4) Detectability Condition (D.C.) 

There exists an operator G G £(5^, if) such that Aq = A + GC 
generates an analytic semigroup Sc(t) = e^ A ^ GC ^ 1 and Sc{t) is expo- 
nentially stable on H, i.e., 

II 5 g(0IIz:(^) < M a e~ WGt , for w G > 0. (16) 

(H-5) In addition to (H-3) and (H-4) there exists a finite-dimensional 
subspace N C W, with dimN < N c such that 

(i) Sp(t)\i C N, for all t > 0, 

(ii) Range GcH, 

(iii) H C V(Ap). 

Moreover, there exist linear maps i : — ► N, n : H — > %l Nc such that 

ni = 7jv c , inx = x for x € N. (17) 

Note that (H-5) implies that ttAfi is a well defined linear map on $t Nc . We 
will show that the system 


w{t) 


n {Ap + GC) i w(t) — irGy(t), 


to (0) — w 0 


( 18 ) 
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u(t) = F iw(t) (19) 

defines a stabilizing compensator for the Cauchy problem (1-2). The follow- 
ing result is a slight extension of Theorem 2.5 in [7] for unbounded inputs 
and outputs. 

Theorem 2.2 If (H-l)-(H-5) are satisfied, then the closed-loop system de- 
fined by (1-2) and (18-19) is exponential stable. 

Proof: Note that without loss of generality we can assume that dim N = N c . 
By Proposition 2.1 it follows that the closed- loop system is a well-posed 
Cauchy problem. Let € W, w a G and z(t), w(t) be defined by (1-2) 
and (18-19), respectively. Since z(t) € W, if x(t) is defined by 

x(t ) = iw(t) — z(t) t > 0, 

then x(t) belongs to W and it is straightforward to show that 

w(t) = x Ajriw(t) + nGCx(t). (20) 

Therefore, 

x(t) = ix Sf(t)iw 0 + j ix Sp(t — s)ixGCx(s)ds — z(t) 

= SF(f)ii«o+ f Sp(t — s) GC x(s) ds — z(t) 

Jo 

= S(t)iw a + f S(t — s)[B F iw(s) + GC x(s)]ds 
Jo 

-S(t)zo+ f S(t - s) Bu(s)ds 
Jo 

= S(t)x(0)+ [ S(t - s)GC x{s)ds, 

Jo 

which implies that x(t) = Sc(t)x( 0). The stability of x(t), w(t ) and z(t) 
follows. 
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3 Optimal Projection Theory 

Consider the steady-state fixed-order dynamic compensator problem, de- 
fined by the infinite-dimensional control system 

z(t) — Az(t) + Bu(t) -f H\7](t) (21) 

with measurements 

y(t) = Cz(t) + H 2 r ] (t). (22) 

The objective is to design a finite-dimensional fixed-order dynamic compen- 
sator 


z c (t) = A c z c (t) + B c y(t) .(23) 

u(<) = C c z c (t) (24) 


which minimizes the steady-state performance criterion 


J(A C , B c ,C c ) *== lim \ f E (i?iz(s) , z(s)) + u(s) T R 2 u(s) 

t-*oo t Jo I 


ds (25) 


where the operators A , B and C satisfy all the assumptions given in the 
previous section and E[ ] is the expectation. In addition, assume that the 
state and measurements are corrupted by a white noise r)(t) in the Hilbert 
space H, with zero-mean Gaussian, Hi € C(H, H), H 2 6 C(H, 3^), Ri £ 
C(H) is self-adjoint and nonnegative definite, and that R 2 is an m x m 
symmetric positive-definite matrix. We assume that the disturbance and 
measurements are independent, i.e., Hi H% = 0, Vi = Hi H* 6 C(H) is 
nonnegative definite and of trace class, and that V 2 = H 2 is 

positive definite. Also, it is assumed that the initial state z(0) = z 0 is 
Gaussian and independent of r/( ). The compensator will be assumed to be 
of fixed, finite order N c (i.e., z c (t) € $t Nc ) and the optimization is performed 
over A c e B c € and C c € 

If one introduces the augmented state space H = H x $t Nc , then the 
closed-loop system becomes a linear system on W. Consequently, define the 
closed-loop operator A : D(A) C H — > H on the dense domain 'D(A) = 
V(A) x by 


A = 


A BC c 
B c C A c 


A 0 

0 0 


+ 


0 BC c 
B c C Ac • 
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Since the operator 


A 0 

0 0 


: V{A) -> H, 


generates an analytic semigroup 


e At 0 

0 /„ c 


t >0, 


then conditions (7)-(8) imply that A is also closed and generates an analytic 
semigroup e M on H (see [10]). To guarantee that J is finite and indepen- 
dent of initial conditions we restrict our attention to the set of admissible 
compensators defined by 

S = {{A c , B c ,C c ) ■ e At is exponentially stable}. (26) 

If {A c , B c , C c ) G S, then there exist a > 1 and (3 > 0 such that 

\\e At \\ < ae-? 1 t > 0. (27) 


Moreover, we know from Theorem 2.2 above that S is non-empty. We now 
state some results found in [1] and [2]. 

Lemma 3.1 If Q and P G C(H) have finite rank and are nonnegative 
definite, then QP is nonnegative semi-simple. Furthermore, if rank ( QP ) = 
N c , then there exist G and T G C(H, %l Nc ) and a positive semi-simple matrix 
M G %t Nc * Nc such that 


QP = G"MT (28) 

TG* = I Nc . (29) 


Proof: Bernstein and Hyland give a complete proof of this result in [1]. Here 
we outline their proof in order to illustrate the form of the factorization of 
QP and to provide a description of the operators G and T. Since Q and 
P have finite rank, there exists a finite dimensional subspace Z C H such 
that QZ C Z, QZ 1 = 0, PZ C Z and PZ L = 0. Hence there exists an 
orthonormal basis for H and in this basis Q and P have the infinite matrix 
representations 


Ql 
0 0 |’ 


P = 


Pi 0 

0 0 ’ 
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where Q\, A € ft rxr and r = dim Z. Consequently, there exists an invert- 
ible $ G ft rxr such that A = is nonnegative and diagonal and 

Qp is nonnegative and semi-simple. If rank (QP) = N c , then it is clear that 
$ can be chosen so that 

A = 



where A G %l Nc x Nc is positive and diagonal. Hence, 


= 0 l] [|] Ml 7 * °] °] . 


QP = 


and if we define G, M and T by 



for any invertible 5 G %t Nc * Nc , then G, T and M provide the desired factor- 
ization and this completes the proof. 


Throughout the paper we will refer to G, T and M satisfying the above 
lemma as a (G - M - V) - factorization of QP. For convenience we define 
£ = BR^ 1 B* and £ = C*Vp 1 C and let I Ne and I H denote respectively 
the N c x N c identity matrix and the identity operator on H , respectively. 
We state Bernstein’s and Hyland’s main theorem which provides a set of 
necessary conditions that characterize the optimal steady-state fixed order 
dynamic compensator for bounded input and output operators (see [1]). 

Theorem 3.2 Let B and C be bounded operators and let N c be given and 
suppose that there exists a controllable and observable dynamic compen- 
sator (A c , B c , C c ) G S of order N c which minimizes J given by (25), then 
there exist nonnegative definite operators Q, P, Q, and P on H such that 
A c , B c , and C c are given by 

A c = T(A-QX-ZP)G * 


(30) 
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B c = FQC*V 2 ~ l (31) 

C c = -R^'B'PV (32) 

for some (G — M — T) - factorization ofQP and such that f with r — G*Y 6 
C(H), the following conditions are satisfied: 

Q : V{A*) -4 V{A) P : V(A) -4 V(A*) 

Q : H -4 V{A) P : H-+V(A*) 

rank ( Q ) — rank ( P ) = rank (C}P) 

and 

0 = (A - tQT,)Q + Q(A - tQT,)* + Vi + tQDQt* (33) 

0 = (A - SFr)*F + P(.4 - EPr) + /?! + r’PSPr (34) 

0 = [(A - DP)Q + Q(A - DPy + QEQ] t* (35) 

o = [(4 - QDyP + P{A - gs) + pep] t. (36) 

Note that these necessary conditions consist of a system of four operator 
equations, including a pair of modified Riccati equations and a pair of mod- 
ified Lyapunov equations which are coupled by the operator r £ C(H ). The 
operator r is idempotent, since r 2 = rr = G*FG*r = G*I nc T = G*F = r. 
In general r is an oblique projection and may not be orthogonal since there 
is no requirement that r be self-adjoint. Moreover, we note that in view of 
Lemma 3.1, Theorem 3.2 applies to (SA C S _1 , SB Cy C c S~ l ) for any invertible 
, since the (G — M — r)-factorization of QP, used to determine 
A Cy B c and C c , is not unique. However, the operator r remains invariant 
over the class of factorizations. An easy computation yields the following 
identities: 

Q = tQ and P = Pr . (37) 

It is helpful to have an alternative form of the optimal projection equa- 
tions to actually compute the optimal fixed-order compensator of the ap- 
proximating finite-dimensional plant. The following result for bounded input 
bounded output operators may be found in [1]. 

Proposition 3.1 IfB andC are bounded, then the optimal projection equa- 
tions (33)-(36) are equivalent , respectively, to 


0 


AQ + QA* + V x - Q^Q + t^QZQt! 


(38) 
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0 = A*P + PA + Ri - PDP + tIPEPtl (39) 

0 = A p Q + QA p + Q^Q-t x QZQt* ± (40) 

0 = A*P + PA q + PEP - t x PDPt x (41) 

where __ 

tj_ = /ii — r, A p = >1 — EP and A 9 = A — QE. (42) 


This form of the optimal projection equations shows that there is a 
connection between Theorem 3.2 and the standard LQG result when dim 
H = N < oo. In this case, we note that the (G — M — T) -factorization of 
QP when N c = N is given by G = T = In and M = QP . Since r = /^r 
and r± = 0, it follows that (38)— (39) reduce to the standard observer and 
regulator Riccati equations. 

To obtain a geometric interpretation of the optimal projection we intro- 
duce the “quasi-full-state” estimate 

m = G*z c (t) e H , (43) 

so that rz(t) = z(t) and z c (t) = Tz(t). Hence, the closed-loop system can 
be written as 


z(t) = Az(t) + BC c rz(t) (44) 

z(t) = T(A + BC c -B c C)Tz(t) + TB c Cz(t) (45) 

where 

B c = QC*Vf 1 and C c = -R^ l B*P. (46) 

This shows that the geometric structure of the quasi-full-order compen- 
sator is dictated by the projection t. Sensor inputs tB c Cz are annihilated 
unless they are contained in TZ(t*) = M(r) L , while rz employed in the con- 
trol input is contained in Consequently, TZ(t) and 7l(r*) are the 

control and observation subspaces of the compensator, respectively. In order 
to modify the previous results so that they will apply directly to unbounded 
B and C operators, care must be exercised to precisely define the weak 
forms of (33)-(36) and (38)-(41). We shall not consider this problem in this 
short note. However, we shall use these systems to guide the approximations 
below. 
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4 Finite Dimensional Approximation 

In general, the optimal projection equations (38)-(41) are infinite dimen- 
sional operator equations. To actually use these equations to compute the 
optimal fixed-finite-order compensator, a finite dimensional approximation 
is needed (see [2] for details). 

Let H n for N = 1,2, be a sequence of finite dimensional linear 
subspaces of H and let V N : H -► H N be the canonical orthogonal 

projections. Let A N G £(H N ), B N G £(» m , H N ), C N € C{H N , 3^), 
R\ G C{H N ) and Vf* G C(H N ) be given and consider the approximating 
system 

z N (t) = A N z N (t) + B N u N (t) + H? V N (t) (47) 

y N (t) = C N z N (t) + H? V N (t). (48) 


The goal is to design a sequence of finite-dimensional dynamic compensators 
of fixed order N c of the form 


i?(t) = A?z?(t) + B?y N (t) (49) 

u N (t) = C c N z N (t), (50) 


which minimizes the performance criterion 


J n (A?, Bf, C t ") = \im \ J‘ e[«^' v ( 8 ) , *»(»)) + , i («) T n 2 „(s)] is. 

(51) 

Now, for each N = 1, 2, • ■ let yfrF j be a basis for H N . Also, for any 

linear operator F N with domain and range in H N , unless otherwise noted, 
we use the same symbol F N for its matrix representation with respect to 
the basis chosen. Let 'I' ;V denote the fc‘ V -square Gram matrix corresponding 

to the basis {<f>? }* =i (e.g., V N - , <!>j) hN \ )• Note that 


(A n Y = (^ Ar ) -1 (A iV ) T 5' iV (B N y = ( B n ) t ^ n (52) 

{c N y = yn N )~ l {c N ) T {t, n ) = b n r^ x {b n ) t ^ n ( 53 ) 

(rf)* = ('I'Wi ) T *" C Z N ) = (54) 



Optimal Fixed’Finite-Dimensional Compensator 


13 


and if we define the k N x k N nonnegative definite matrices 


Qg 

def 


pN def qNpN 

Qo 

def 

QW 1 

pN def yNpN 

v 0 N 

def 

vfw")- 1 

< d i' 


def 

B N R^ l (B N ) T 

(C" 


N ' T vr l c s 


then the matrix equivalence of the operator equation (38)-(41) become 


0 = 

A n QH + QZ(A N ) T + V 0 N - Qo'BqQq 



+ t jl Q o^o Qo ( T i r ) T 

(55) 

0 = 

(A N ) T P 0 N + P 0 n A n + <- Pq* Eq Pq 



+(t?) T Po N Z»P 0 n t? 

(56) 

0 = 

+ <? 0 "(<) T + Qo^oQo 



~ T ± Qo^o Qo ( r ] v ) T 

(57) 

0 = 

(<) T ^ + + PfzgPp 



~( T ± ) T P^ ^o Pq 1 t l • 

(58) 

The approximating 

optimal dynamic compensator (A^B^C^) of order 

N c is then given by 



A? 

= r*{A N - - Eq Pq*)(Gq) t 

(59) 

B? 

= Ti ?Q$(C n ) t V 2 - 1 

(60) 

c? 

= -i^ 1 (B Ar ) T F 0 Ar (r^) T 

(61) 


where Tff, Gq €jR N ‘ xkN and G provide a (Gq - M? - rjf)- 

factorization of Qq Pq. 

We turn now to an example. Consider Burgers’ equation, with Neumann 
boundary control given by 


d , x 

z(0,x) 


d_ 

dx 


z(t,0) 


Sr d 

e-^z(t,x) - z(t, x )— z (t,x), 0 < x < 1, <>0 (62) 
z 0 (x) (63) 




(64) 


^2(t,l) =U 2 (<), 
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and observations 


Vi (0 = z(t,0) (65) 

1/2 «) = z(t, 1), (66) 


where e = rr- > 0 and Re is the Reynolds number. Initially, we consider 
the linearizea Neumann boundary control problem 

z(t,x ) = ej^jz(t,x), 0 < x < 1, <>0 (67) 

z(0, x) = zo(x) (68) 

J^z(i,0) = -u x (t), -^-z(t,l) = u 2 (t). (69) 

We will apply the linearized feedback control laws constructed from this 
model to the nonlinear Burgers’ equation. System (67)-(69) can be placed 
into the standard state space framework by defining the operator A t on 
H = L 2 ( 0,1) by 

A t 4> = e<f>" (70) 

for all <f> € V(A ( ) = {<j> € H 2 { 0, 1) : ^(0) = ^(1) = 0}. Define W. = V* = 

H?( 0, 1) = D(A^) and let B : Si 2 V be defined by B = AAf where A = 

—A t +wl and we assume that w is not an eigenvalue of A t with homogeneous 

Neumann boundary conditions, so that A is boundedly invertible on L 2 (0, 1). 
The Neumann map M is defined by the boundary system given in [11, pages 
53-56]. Let C : W -> 3ft 2 defined by 


C<t> = 


m 

m 


(71) 


The boundary control problem (67)- (69) can be represented by a differential 
equation 

4- = A ( z(t) + Bu , x(0) = z 0 (72) 

at 

y(t) = Cz(t). (73) 

It is well known that A t generates an analytic semigroup S(t) on H. More- 
over, the spectrum o(A<) of A t consists of all eigenvalues A„, n = 0, 1, 2, • • ■ 
given by A„ = -en 2 n 2 and for each eigenvalue A„ the corresponding eigen- 
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function <f>„ is given by 

<i>o(x) = 1 <j>n(x) = y/2cos(nnx), 0 < x < 1. (74) 

One can easily verify conditions (7)-(8), by taking 7 = £ (see [11]). It is 
straightforward to show that (H-l)-(H-5) are valid. 


5 Numerical Results 


Now, we formulate a specific approximation scheme for the boundary control 
problem (74). For each N = 2, 3, • • • let divide the unit interval [0,1] into N 
equal subinterval [xi,x<+i], x, = i = 1, 2, • • • , N + 1. Let H N = Span 

KLo where are the standard hat functions defining continuous 

piecewise linear splines (see [3]). Note that k N = dim H N = N + 2 and let 
the approximate solution z N (t,x) of z(t,x) for equation (72)-(73) be given 
by 

N + 1 

z N (t,x)= £ Z?(t)h?(x) (75) 

1=0 


for some zf* (t) € 3?, i = 0, 1, • • • , N + 1. Standard finite element approxima- 
tions yield the ODE system 


= A"z N {t) + B N u{t), z N (0) = zg (76) 

y N (t ) = C N z N (t) (77) 

where the matrices A^' , B N , C N can be easely computed by using the 
Ritz-Galerkin approximation. 

For our numerical example, we set a = ^, the initial condition z Q (x) = 
sin(irx), r\ = v\ = 1 and r 2 = v 2 = 10 -3 . Also, R\ = r\I H , R 2 = r 2 / m , 
V\ = v\Ih, and V 2 = v 2 /(. Therefore, it follows from Section 3 that R° = 
ri^! N and = iq^^) -1 where is the Gram matrix. In this numerical 
example we will compare the approximating optimal LQG (i.e., N c = N + 2) 
with the dynamic compensators of various order N c . The optimal projection 
equations (55)-(58) were solved using the homotopic continuation algorithm 
described in [16]. The approximating controllers defined by the linear fixed- 
order compensator (B? and C?) were applied to Burgers’ equation (62)- 
( 66 ). 
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We note that N = 32 produces converged optimal LQG designs. Hence, 
the reduced order compensators were tested on both the linear and nonlinear 
problem using the N = 32 order finite element model. 

In the full order case N = 32 and N c = 34, the converged feedback and 
observer functional gains are given in Figures 1 and 2, respectively. Since 
we are controlling the flux at each end point x = 0 and x = 1 , we have 
two feedback functional gains, the one plotted with solid line is the flux 
control gain at the origin and the one plotted with dashed line is the flux 
control at the end point x = 1. Similarly, since we are sensing the flow 
at the origin and the end point, we have two observer gains (solid line for 
observer gain at the origin and dashed line for the observer gain at x = 1). 
Next, we applied the full order controller to Burgers’ equation resulting in 
the nonlinear closed-loop trajectory given in Figure 3. 

In the fixed-order case, we considered the accuracy of the impulse and 
step responses of the various reduced order compensator designs compared 
to the corresponding responses of the full order LQG design. Figure 4 il- 
lustrates the linear closed-loop impulse response for the full-order LQG and 
reduced order compensator (of order N c = 16) designs. The impulse re- 
sponse of the linear closed-loop system for the 16th-order compensator is 
in perfect agreement with the LQG response. Note that in Figure 4 we see 
only one plot for both designs because both plots are essentially the same. 
Similar trends are seen (Figure 5) in the comparisons of the step responses 
(for the same design case) with the corresponding LQG responses. 

For the nonlinear closed-loop response, the 16th-order compensator was 
applied to Burgers’ equation and we see (in Figure 6) excellent agreement 
with the full order closed-loop trajectory response. Hence, replacing the 
32nd-order optimal LQG controller by a 16th-order compensator produces 
a closed-loop system with minor performance degradation. 

We also compared the performances of the closed-loop system of the 4th- 
order compensator with the full order LQG responses. Figures 7 and 8 are 
the impulse and step responses of the linear closed-loop system, respectively. 
If one compares these responses with the corresponding responses for the full 
order LQG controller shown in Figures 4 and 5, then it is clear that the 4th- 
order compensator performs almost as well as the full order LQG controller. 
Similar comments hold for the nonlinear closed-loop responses. For example, 
the 4th-order compensator response (Figure 9, solid line) compares well to 
the LQG response (Figure 9, Dashed lines), especially after time T = 1.0. 
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Figure 9: Closed-Loop for the LQG & Fixed-Order 
Cases, N=32, N c = 4 


6 Conclusion 

The purpose of this note was to show that finite dimensional dynamic com- 
pensators could be used to control a nonlinear partial differential equation 
without significant loss in performance. Although there is considerable the- 
oretical and numerical work for bounded input and bounded output oper- 
ators, numerical results for the unbounded control and observation opera- 
tors are not as fully developed as the theory. For example, several authors 
have considered questions of existence of stabilizing dynamic compensators 
(even for nonlinear plants [6]) for boundary control problems. However, 
approaches, such as the optimal projection method, that result in a “com- 
putable” fixed order compensator have not been applied to more general 
boundary control problems. Although the numerical results presented here 
show that the optimal projection method can produce excellent designs for 
problems with boundary control and observation, there are a number of the- 
oretical and numerical issues that need to be resolved in order to extend this 
approach to practical problems of this type. 
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